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ABSTRACT 

By observing the transit of various cloud features across the 
Jovian disk, Terri le and Westphal (1977) have constructed lint) -darkening 
curves for three regions in the 4.6 to 5.1 um band. Several models cur- 
rently employed in describing the radiative or dynamical properties of 
planetary atmospheres are here examined to understand their implications 
for 1 ini) -darkening. The statistical problem of fitting these models to 
the observed data is reviewd and methods for applying multiple regression 
analysis are discussed. Analysis of variance techniques are introduced to 
test the viability of a given physical process as a cause of the observed 
limb-darkening . 

The intermediate flux region of the North Equatorial Belt appears 
to be in only modest departure from radiative equilibrium The limb- 
darkening curve for the South Temperate Belt is rich in structure and 
cannot be satisfactorily ascribed to any single physical mechanism; a 
combination of several, as yet unidentified, processes is likely involved. 
The hottest areas of the North and South Equatorial Belts exhibit limb- 
darkening curves that are typical of atmospheres in convective equilibrium. 
In this case, we derive a measure of the departure of the lapse rate from 
the dry adiabatic value (n~1.68), which furnishes strong evidence for a 
phase transition at unit optical depth in the NEB and SEB. Although the 
system NH^ cannot be entirely ruled out, the freezing of an aqueous 
anmonia solution is shown to be consistent with the parameter fit and 
solar abundance data, while being in close agreement with Lewis' (1969a) 
cloud models. 


I. INTRODUCTION 


By applying a combination of radiative transfer and statistical 
techniques to infrared observations of Jupiter, we can enhance our 
understanding of the radiative and dynamical processes that control the 
make-up of the Jovian clouds and the deeper atmosphere. Since Jupiter's 
atmosphere is widely believed to be significantly stratified in its 
spectroscopically active components, observations made over limited 
wavelength regions permit us to look down at levels in the atmosphere 
where those components play an active role. Moreover, by analyzing the 
intensity observed in different regions of the disk, some of the under- 
lying physics of Jupiter's belts and zones can be revealed. 

The 5 ym region is of particular interest since it is transparent 
to the abundant Jovian absorbers, gaseous hydrogen, methane and anmonia. 
Gillett et al . , (1969) observed that the 5 ym brightness temperature , 
averaged over a large part of the Jovian disk, was approximately 230°K. 

This value is much higher than most workers had expected and, since it 
corresponds to deep atmospheric levels, was a stimulus to further in- 
vestigations. Westphal (1969), in observations of the North Equatorial 
Belt, showed that the 5 ym flux was coming from localized hot spots with 
brightness te m p er atures >300°K. He concluded that, if it is assured that 
the cloud layer was near the top of the convective zone, the radiation 
was coming from below the clouds. Keay et al. , (1973) and Westphal et al. , 
(1974) produced high resolution maps of Jupiter confirming the existence 
of localized hot spots. In addition, they observed a correlation of 5 ym 
features with visual features in color photography. High thermal flux 
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seemed to con® from "blue” or "purple" regions whereas "orange" or "red" 
regions were not sources of intense 5 m radiation. They speculated 
that this dichotomy was likely due to the absence or presence of middle 
altitude red clouds. Sagan (1971) had earlier proposed that the blue 
coloration vas due to Rayleigh scattering at roughly the 1 bar level 
when our view is not impeded by intervening clouds of red chromophores . 
Westphal et al , concluded (by studying the flux emanating from the shadow 
of lo as it passed across the Jovian disk) that the 5 ym flux was not re- 
flected or scattered sunlight; but a genuine feature of radiative sources 
deep within Jupiter's atmosphere. 

Tb better understand the nature of sen® of Jupiter's 5 ym features, 
Terrile and Westphal (1977) measured limb-darkening by selecting a region 
of interest and measuring its brightness as it rotated around the planet. 

In particular, they observed the hottest emitting areas of the North and 
South Equatorial Belts (with brightness temperatures of about 250 to 255°K) , 
several bright 5 ym areas in the South Temper* t? Belt (with similar 
brightness tenperatures) . and intermediate flux regions in the North 
Equatorial Belt (with brightness tenperatures near 240°K) . This mtailed 
the use of many different images of Jupiter taken over several hours in 
order to construct one limb-darkening curve. Notably, this technique did 
not suffer from the smearing effects of longitudinal inhomogeneities that 
would result from generating Hub-darkening curves from thermal maps. 

In parallel with these infrared observations, the atmosphere and 
conposition of Jipiter was undergoing extensive study. Lewis (1969b) es- 
tablished that if the Jovian atmosphere possessed the solar abundance of 
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water and amnonia, the clouds were dominated by an aqueous anmonia solu- 
tion while the topmost cloud layer was solid anmonia. He also showed 
that If sulfur were present in solar abundance , NH^SH would form an im- 
portant cloud layer. The infrared properties of liquid and solid watt- 
(Irvine and Pollack, 1968; Robertson and Williams, 1971) and arronia 
(Robertson mid Williams, 1973; Robert sen et al. , 1975) have been investi- 
gated at 5 \m and are known to have very large absorption coefficients. 

These 5 van observations, together with predictions obtained from 
models of Jupiter ' s atmosphere and the infrared properties of its con- 
jectured constituents, provide a compelling reason for analyzing the ob- 
served limb-darkening of various regions of the planet. In this communica- 
tion, we present a combined radiative , dynamical and chemical model that 
reproduces the observed limb-darkening curves. 

II . RADIATIVE AND DYNAMICAL 1CDELS 

The plane-parallel approximation to the equation of radiative 
transfer is 

v ft I (T,u) - I (t,p) - S(t) (1) 

where arccos w is the angle between the line of sight and the local plane- 
tary normal, t is the optical depth, I is the intensity of the radiation 
and S is the source function. The optical depth is defined, in differen- 
tial form, to be di = - <dz, where k is the extinction coefficient, and 
z the altitude. All quantities in (1) are considered to have been modu- 
lated by the spectral response function of the InSb detector employed by 
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Turtle and Vtestphal and integrated over the instrument ’ s 4,6 to 5.1 van 
bandwidth. Formally, (1) can be integrated to give the limb-<.iar kening 
fine t ion. 

1(0, u) - J S <t) if 1 exp (-t/u) di (2) 

(Chandrasekhar , 1960). This equation is ideally suited to our analysis 
since Terrile and Vtestphal (197?) have evaluated 1(0, u). Thus, by in- 
verting (2), we can detemdne the source function S(0- 

In general, the inversion is unique only if the limb-darkening 
funct ion 1(0, y) has a known functional form. In the case of discrete 
data (particularly data contaminated by noise) , the inversion is not 
unique and we trust choose one of two approaches. In one approach, we cal- 
culate an approximate inhesion kernel, K (t,u), such that 

3(0 * ^ K(ui) l(0,u) du . (3) 

(See the recent review article by Parker, 1977, for a description of 
’’generalised inverse theory."), Gooidning (2) and (3). we require that 

AO,)’) kci.uiu' 1 exp (-iVt'.)di. <A) 

".proximate” the Dirac delta function, 6 ( \ - \ ' ) ; that Is, that the integral 
of over \ or \ ’ is unity and A(t . i * ) is strongly peaked when t’ ap- 

proaches t . The inversion is ;aade complete hv specifying the approximate 
source function and an os t Invite of the width of A(i , i ’) for different values 
of i (which, in turn, provides a measure of tl <e characteristic convolutional 
smoothing evident in tlx approximate source function). Orton (1977) re- 
cent Ay aiployed this technique in recowritv, tlx mean Jovian tenperat\a*e 
structure from spectrally resolved thortm 1 radiance hit a . This method, 
txwevor, provides no direct insight into the physical processes that are 
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the source of the infrared radiation and, mor e o ver , has several mathe- 
matical deficiencies that are apparently not well-known (see Appendix I) . 

A less general but more physically motivated approach is to con- 
struct several radiative and dynanical models and obtain their corres- 
ponding source and limb-darkening functions. These models will depend, 
often nonlinear ly, on a small mnfcer of parameters. By employing nultiple 
regression methods , we then obtain miner i cal estimates of the parameters 
that are in a statistical sense most likely. By then employing analysis 
of variance techniques, we can assess \&ether the residual errors in the 
model fits are compatible with the experimental noise. This approach can 
demonstrate directly that a given physical model could be responsible 
for the observed lint)- darkening while other models mist be rejected. It 
is important to note that the fit obtained is most accurate for t~l, the 
vicinity of the cloud tops and the region of greatest physical interest. 

The reason for this is clear from (2) . The source function for t>0 depends 
strongly on 1(0, y) measured near the lint) where instrumental accuracy is 
least. The source function for t»1 is strongly attenuated and could vary 
significantly without seriously affecting the observed limb-darkening 
function. Because of the intuitive value of the method, we shall confine 
our attention to this technique and turn now to a discussion of models. 

A. Power Series Expansion Model 

Although the power series expansion 

1(0,0 -X «n l ' n < 5 > 

has no direct physical interpretation, we choose to include it for several 
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reasons. The truncated expansion was useful in an analysis of the limb- 
darkening of Venus (Goody, 1965; Newman, 1975), It can provide an es- 
timate for the scatter in the data due to noise that is required in the 
analysis of variance. Finally, assuning that l(0,u) is analytic and 
regular, the inversion of (2) can be performed directly giving (5) where 
the source function may be written 

s <’> -Xv n/n: »> 

In one instance, the well-known Eddington approximation, the truncated 
power series expansion is of special interest; 

1(0, p) - 1(0,0) (1 4- 3A0 . (7) 

This limb-darkening function results if the flux over the corresponding 
frequency passband is conserved . If the flux over the entire frequency 
spec trim is conserved, we have radiative equilibrium. 

B. Convective Equilibrium Model 

In our terrestrial experience, clouds are very efficient infrared 
absorbers of solar radiation as well as heat from the surface (additionally, 
in Jupiter's case, heat generated internally), A warmed parcel of gas 
will rise and adjust its pressure to that of its surroundings , the pressure 
of which varies according to the equation of hydrostatic equilibrium, 

dP/dz - -eg (8) 

where P is the pressure, e the atmospheric mass density, and g the local 
gravitational acceleration. Thermal conduction times are very slow com- 
pared with dynamical times and nay be neglected. Thus, the resulting 



behavior of the parcel of gas is adiabatic and the pressure in the 
parcel behaves according to 

P “ P Y (9) 

vrtiere p is the parcel’s mass density and y is the ratio of specific 
heats. 

Assuring that the absorbers responsible for the extinction are 
well-mixed with the principal atmospheric constituents, we can express 
the extinction coefficient as 

k * oxp/m^ (10) 

*&ere a is the cross-section to absorption, x is the mixing ratiu of the 
absorbers to the principal atmospheric constituents , and m^ is the mean 
mass of a single absorber molecule. Combining (2), (8), and (10), we ob- 
tain a linear dependence of pressure on optical depth, namely, 

dP/dT = m^hx. . (U) 

By terrestrial analogy, we expect clouds to have a fairly sharp top at a 
level characterized by a temperature and a ores sure P . The discon- 
tinuous boundary may result from the transition from convective to radia- 
tive equilibrium or from phase changes. Defining t to be zero above the 
cloud top (where we assume there to be no significant absorption) , 
equation (11) is integrated as 



Pollack and Sagan (1965) derived a similar expression for the Venus atmos- 
phere. In their case, however, the absorber was the principal atmospheric 
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constituent and x-1. 

From (9) and the ideal gas law, we find 



(13) 


Let us now assure that the clouds radiate as a black-body. By inte- 
grating the Planck function from 4,6 to 5,1 \m (assuring a relatively 
uniform response in the InSb Interference filter enployed by Terrile 
and Westphal) , we obtain the approximate power -law dependence 


B 5^ • < 14 > 

Typical values for the exponent n are 12.3 and 11.6 for temperatures 
of 240°K and 250°K, respectively. Confining (12) , (13) , and (14) yields 


the intensity dependence on optical depth 


B 5um (T) ’ B 5mn (0) 


1 + 




1 


n Cr-i? 

Y 


(15) 


If we combine the equation of state with (8) and (9) , we obtain 


the dry adiabatic lapse rate 


where 

atm 

constant. 


dr 

ai 


m 


atm y-1 


~if T * (16) 

is the mass of an atmospheric constituent and k is the Eoltanann 
(For an inhomogeneous atmosphere, the meaning of an "atmospheric 


constituent" may be ambiguous. We define such a pseudo-particle as being 
characterized by a mnber-density weighted average of each component . In 
the terrestrial case, the mass of such a fictitious constituent wuld then 
be 0.78 the mass of 0.21. times the mass of O 2 and 0.01 the mass of a 



trace constituent . ) This equation indicates how the temperature of 
our parcel of gas decreases as it rises. 

Suppose that the parcel of gas has a trace of a substance that is 
undergoing a phase transition, for exarple water vapor in the terrestrial 
at mosphe r e . As the vapor condenses, it evolves heat and precipitates 
out of the parcel. Because of latent heat it deposits in the parcel, 
the lapse rate is reduced. Li particular (see Hess, 1959) the term 
(y~1)/y in (15) and (16) should be replaced by (y-1)/hy, whore n is given 


by 


1 + 


el/ 

5L 


w 

V 


1 + 


Lw 

ET 


(17) 


and Cp is the mean heat capacity of an atmospheric constituent, w the 
mixing radio of condensates to atmospheric constituents, r the ratio of 
the molecular weights of the condensates to that of the atmospheric con- 
stituents, and L the latent heat evolved by a single condensing moleoile. 
Equation (16), when corrected for the condensate, defines the wet adia- 
batic lapse rate. Similarly, U5) becomes 

B 5w/> ; B 5/ 0) 

and describes the effective black body intensity as a function of optical 
depth for a wet adiabat . 

Let us assure that the source function S(i ) may be appro'dmated by 
the effective black body intensity. For simplicity we write (18) as 

S ( r ) - a(l+bi ) C 


1 + 


^abs 

tbdV 


n y 


(18) 


(19) 
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where a, b, and c correspond to appropriate terms in the previous 
equation. Then, employing (2), we find 


1(0, y) • a(by) a exp 





where r is the incomplete Gwnna Function (Abranrxrttz and Stegun, 
1965) . A convenient formula for evaluating (20) is 


1(0, y) - a(by) c exp 


(by) 


-1 


{r < 1+c >-n£o Sftd 5 


-1 


n 


where r here denotes die complete Garnna Function. 


( 20 ) 


( 21 ) 


C. Cloud and Intermediate Zone Models 


These models were enployed by Terri le and Westphal (1977) in analy- 
zing their data. Their cloud model describes radiation from an optically 
thick, hot cloud deck passing tlrrough an optically thin, warm, emitting 
layer. This may be represented by 


1(0, y) 


exp(-T w /y) + \ 


■ 

l-exp(-T w /u) 


( 22 ) 


where and 

warm layers respectively, and t is the optical depth of the warn layer. 
Radiation from an ootically thick intermediate cloud deck passing through 
a cold absorbiiv; layer (I , intermediate model) may be represented as 


are the black body radiation emitted from he hot and 


1(0, y) * Bj exp( -Tt/u) 

Both of these models can be expressed in the form 


(23) 


1(0, u) - a + b expC-x'/w) 


(24) 
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where the corresponding source finction is given by 

oW 

,'5, • C 

We do not expect these models to give particularly good fits because 
they require that each cloud and absorbing layer have a uniform tem- 
perature distribution. They are included, however, because they pro- 
vide some insight into the nature of lint -darkening functions produced 
by atrospheres with relatively little temperature structure. 



D. Thin Shell Model 

It is also of interest to examine the limb-darkening function of 
an atmosphere characterized by the opposite extreme: an extremely hot, 

very thin emitter embedded at t q in a warm, absorbing atmosphere. The 
source function used is 

S(t) = a +b6(x-x o ) (26) 

where 6 is the Dirac delta function. The corresponding linb-darkening 
function is 

1(0, y) = a + by" 1 exp(-T Q /y) . (27) 

Although this model has no known physical counterpart;, it provides a 
useful measure of the impact of strong temperature variation and pro- 
nounced thermal structure on the limb- darkening function. 

Although our list of models is small, it describes a wide spectrum 
of behavior. The power series expansion describes virtually any con- 
tinuous, smoothly varying limb-darkening function. Moreover, it can be 



used to provide an estimate of the scatter In the data due to noise. 

A special case of the polynomial is the Eddington approximation which 
provides an accurate representation of the limi>- darkening curve when 
the radiative flux is conserved over a given pasaband, a hint of 
possible radiative equilibrium. We expect radiative equilibrium to be 
a dominant feature of high-elevation clouds (not under the influence 
of direct heating from the planet’s interior) if dynamical effects are 
unimportant (a situation found to be the cas. on Venus by Newman, 1975). 

CP the other hand, we expect that deeper clouds are dominated by dy- 
namical effects and, from the suspected caiposition and temperature 
range of these clouds (Lewis, 1969b), subject to phase transitions in 
their spectroscopically active components. For this reason, we haw em- 
ph zed the derivation of the convective equilibrium model. Finally, 
as measures of the degree of therr.nl structure for the observed limb- 
darkening, we also consider Terrile and Westphal's (1977) cloud and inter- 
mediate zone models as well as a thin, extremely hot shell model. The 
source and liirb-darkening functions for the different models are tabulated 


below. 
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Model 

Source Function 

Limb-Darkening Function 

Pcwer- Series 
Expansion 

Jo 

Jo*n“ n 

Convective 

Equilibrium 

a(l-Rn) c 

a(bp) a exp (bp)” 1 P 

* 

1+c, (bp) 1 

Cloud and 

Intermediate 

Zone 

a, 0 -tV 

a + b, t'— t 

a4b exp(-i ' /p) 

Thin Shell 

a+b5(T-i o ) 

abb p” 1 exp(-r o /p) 


Table 1. Lurb-Daritening Propet ties of Various Models 

III. STATISTICAL METHODS 1 AND ANALYSIS 
Suppose we have N measurements, at various zenith angles (arccos 

p^, i = 1 N) of the limb-darkening function (which we denote by 1^) . 

For simple cases, a model may be considered composed of a linear combination 
of M different functions of p, say f. (p) ; j = 1 M. For example, a 

1/ No single reference provides an adequate survey of this problem. 

Ralston (1965) reviews some of the nimerical problems associated with 
least-squares techniques. Jenkins and Watts (1968) examine the theory of 
maximum likelihood estimators and Gaussian least squares as well as pro- 
viding some remarks on nonlinear problems. Graybill (1968) considers 
the general linear model and some statistical tests of confidence. 
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? M-l 

power aeries expansion employs the finctions 1, p, p , u 

and we make the identification f^(u) * p^“ . Nonlinear models will be 

treated later in this section. In addition, we assure that there is an 

additive Gaussian error noise component, r, , whose mean vanishes and has 
2 

a variance of o . Therefore, v* write 


M 

I t - a j fjO^) + e v i - l f . . . . , N (28) 

vrtrere the are linear corbination coefficients. We consider the error 
to be Gaussian distributed, a reasonable assumption from the Central 
Limit Theorem. Systematic errors, notably those due to calibration, are 
not Gaussian, and are often intractable. 

Since the errors r ^ defined by (28) are Gaussian distributed, the 
probability associated with the estimates of the a^ coefficients varies as 


exp 


N 


i”i 


4 _ f 4 . 

r i /l ' 


Therefore , the most probable choice of the coefficients is that which 
minimizes 


N 

u - >: 

i=l 

The process of finding the values for the a ^ coeff icient.s is called Lite 
maxim m likelihood method ark! is equivalent to the method of least 
squares. In this, the linear case, it is also known as the "multiple 
regression" nixiel . We therefore require that the derivative of U wit h 
respect to each aj coeff ieient vanish, yielding the normal equations 


h - 


j \ -jW 


/ 
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N 

r 

i-1 




M 

V 

k-1 


\\ 


N 

v 

i-1 


f j (M t ) 




(30) 


V i 

This set of linear equations is characterized by a matrix whose j 
component is given by 


<f j' f k ) 7 j L f j (,1 i ) f k (li i> • < 31) 

The matrix is symmetric , send. -positive definite and the system of 
equations has a unique solution unless one of the functions f^ (u^) , 

evaluated at each i-1 N, could be represented by a linear 

combination of the remaining functions and is thus redundant. 

Although superficially sinple to solve, the system of linear 
equations (30) is numerically ill-conditioned. For example, fin eighth 
degree polynomial fit to an arbitrarily large data set will result in 
the loss of twelve significant places of accuracy if a direct method 
(e . g . . Gaussian elimination with pivoting) is usedl To reduce this 
source of computational error, one should use Gram-Schnidt orthogpnaliza- 

t ion of the functions fj (u) , j - 1 M, with respect to die inner 

product operator defined in (31). 'Hie resulting matrix, associated with 
the normal equations will then be diagonal fu id the system's solution will 
then be trivial to obtain. A relatively recent innovation in solving 
least-squares problems is the technique of singular value decomposition. 
Although functionally equivalent to the Gram-Sclmidt proctxiure, it is 
somewhat faster in execution. Moreover, unless pt-sed in a certain form, 
the Gram-Schnidt procedure is susceptible to numerical instabilities. 
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Singular value decomposition is a very complex procedure but is de- 
scribed in detail in Lawson and Hanson (1974) and Forsythe et al . , 
(1977). Also, both texts contain tested At SI Standard Fortran pro- 
grams. 

Although the errors associated with an ideal experiment are in- 
dependent of each other, the errors estimated by (28), where the a^ 
coefficients satisfy (30), are not. In fact, combining (28) and (30), 
we find 


N 

£ e.f,(u.) - 0 for j - 1, . . . ,M . (32) 

i-l 1 J 1 

So, although there are N values of e^, equations (32) introduce M con- 
ditions or constraints and we are left with N-M degrees of freedom. 
Moreover , if we calculate the expectation value of U (the sum of the 
residual variances) defined by (29), we can shew that 

N 2 2 

<U> - E <e • > - (N-M) 0 ^ (33) 

i-l 1 


The least squares estimation processes introduce a small bias (which 
vanishes as the rairber of data points becomes arbitrarily large) due to 
(32) into our estimates of the fu . By increasing the lumber of a., co- 
efficients, equations (32) shew that we reduce the noise level in each 
s ^ until the number of coefficients M equals the lumber of data noints N 
and all vanish (i.e., the fit is exact) . This, qualitatively, is the 
result show bv (33) . can aprroxiimte the latter by writing 


Jl - 1 
0 _ tra 


N 

v r * 

i*l E 1 


(34) 


(This result is exact only if we replace by <e^ >.) Suppose, for 
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example, that we knew our liirb-darkening function has an exact 

representation, apart from noise, as a polynomial of degree H* , 

2 

Using (34) , we can estimate o , If we fit the data with a poly- 

2 

noraial of degree h>H’ , our estimate of o from (34) will retain ap- 
proximately the same because the decrease in the noise level is 
exactly compensated for by the ctencsninator N-M. However, if we use 
an estinate of the polynomial degree M<M' , we will find that the 


residuals contain not only noise information but liirb -darkening 

2 

information as well, and our estimate of o will be too large. 

2 

Using (34) we can new define a x variable with N-M degrees of 
freedom, namely 


2 

\ 


N 

y 

i“l 



(36) 


2 

If o is knewn (i.e. , we have an absolute estimate of our sources of 

error) , we can enploy the usual con f i dene e - 1 eve 1 tests . 

2 2 

If a" is not knewn, the x test cannot be used. Ihe problem of 

finding the polynomial degree M* is then complicated by the fact that 

2 

(34) is an approximation. In practice, we find that o decreases with 

9 

increasing M, until M equals M' . For larger M, the estimate of o 
tends to oscillate aroirtd a constant, making the task of identifying M' 
very difficult. 

By msxlifying a technique developed by Akaike (1969) in application 
to autoregressive decomposition, we can construct a variable that will 
better equip us to estimate M' . We note that, as M increases , the a^ 
coefficients adapt to the properties of the noise in that experiment 
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until the ranker of coefficients and data points are die same and no 
residual errors regain, let us suppose that, next door to Terrile 
and Westphal , there was a carpeting group using equivalent instrumenta- 
tion making die same observations at the same zenith angles. The 
underlying link-darkening function would be the same for both groups 
but the noise observed presumably would not. (The errors of course, 
would be drawn from the same statistical population.) We then ask how 
well the Sj coefficients conputed for Terrile and Westphal’ s observations 
would match their rivals’ data. That is, if their conpetitors observed 
intensities i-1, . . . ,N, how large would U A be, where we define 


U, 


N 

Z 

i-1 


, M 

X i ' S jW 


A straightforward but tedious calculation reveals that 


(36) 


<U A > - (two O 2 -^ <U> (37) 

when M>M* , As M increases beyond M* , (as an approximation to < U A >) 
increases because the surfeit of coefficients are adding to sane of 
the noise level observed in the first experiment. This variable there- 
fore provides a test of the universality of the fit. In practice, we can 
only estiimte < v by evaluating £(N+M)/ (N-M)J U, for increasing values 
of M. The resulting locus of points is parabolic in character with M’ 
corresponding to the mininun. Since we have obtained 13^ (and not 
the points may oscillate, but the uncertainty in M’ is characteristically 
reduced. Finally, knowing M’ , we can estimate o from (34). This value, 



19 


however. Is not accurate enough to permit anything but the crudest \ 
test of significance. 

As an illustration of these methods, we tabulate seme relevant 
quantities for polynomials fitted to Terrile and Vtestphal's observa- 
tions. 


Polynomial 

Case A 

Case B 

Case C 

Equatorial Belt 

South Torperate 

Intermediate 

Results 

list Spots 

Belt 

Flux (NEB) 

N 

100 

50 

24 

U 

1st degree 

0.1754 

0,2181 

0.04681 

2nd degree 

0.1370 

0.1343 

0,04677 

3rd degree 

0.1358 

0.1327 

0.03832 

4th degree 

0.1323 

0.1207 

0.03667 

5th degree 

0.1316 

0.1097 

0.03586 

<5 

1st degree 

0.04230 

0.06741 

0.04613 

2nd degree 

0.03758 

0.05346 

0.04719 

3rd degree 

0.03761 

0.05371 

0.04377 

4th degree 

0.03732 

0.05179 

0.04391 

5th degree 

0.03742 

0.04994 

0.04463 

U A 

1st degree 

0.1825 

0.2363 

0.05532 

2nd degree 

0.1455 

0.1514 

0.06010 

3rd degree 

0.1471 

0.1558 

0.05365 

4th degree 

0.1462 

0,1475 

0.05598 

5th degree 

0,1484 

0,1397 

0.05976 


Table II. Fewer Series Expansion Parameters 


In Case A, the Aka ike criterion would clearly select a quadratic fit. 
The oscillation in U A is not a hindrance here. In Case B, the mtnlnun is 
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readied for a fifth degree polynomial. However, that fit (as well as 
the 7^ and 8^ order fits) extrapolates to a negative intensity 
at the linb and must be disregarded. We consider, accordingly, a quart ic 
polynomial to be appropriate. Case C is same*bat anbiguous because of 
the strong oscillation in U A< Although M-3 is a minimum, the corres- 
ponding coefficients would provide a source function that was negative 
at t ■ 1.41 and mist be excluded. Hence, the first degree polynomial is 
selected. It is important to note that, in this case, physical and not 
statistical considerations resolved die degree of the polynomial fit. 

Honlinear models are significantly more difficult to fit and 
analyze than their linear counterparts. Instead of linear combination 

coefficients aj . « will <*nloy parsers «, . J - 1.2 M, so that we 

can parallel equation (28) by writing 

1^ *■ F ; a ^ , . . . , a^) + e i (38) 

where the function F describes our model (such as the convective or cloud 
models) . We define the residual variance U by 

U 

We perform a variation of the parameters so as to minimize U, in com- 
pliance with the maxiimn likelihood principle. Unlike the linear case, 
there may be several minima and a global search must be performed. 

The maximum likelihood estimates of the parameters satisfy the 
normal equations (derived by differentiating U with respect to a^) , 
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N 


3F 


N 


± l x h HJ (y i ; *i""> " t ; x F(u i ; *v 


Hj; (p i ; a i< 


j • i h 



( 40 ) 


The minimization of (39) or, alternatively, the solution of (40) is a 
very difficult computational problem (consider the convective model (21) , 
for example) , A survey of this problem may be found hi Luenberger (1973) , 
As in die linear case, we have N measures of the error (38) m id M 

constraints, equations (40), By linearising F(u; ^ a^), we can 

demonstrate the approximate validity of (33) -(35) in the nonlinear problem. 
Thus, once we have obtained the roaxinun likelihood estimate of the 
paraneters, the statistical method of analysis is much the sane as before. 
Because of nonlinearity, there is no direct analogue to Aka ike's criterion. 
In comparing the residual variances U for different models, w re- 
quire, following (34), that they have the same nurber of parameters . thus, 
the intermediate zone model may be compared with a first degree polynomial , 
and the convective equilibrium, thin shell or cloud models with a second 
degree polynomial. Higher order power series trust be treated on an indivi- 
dual basis. 

The residual variances U of the models considered are tabulated 


below. 


Variance 

Case A 

Case B 

Case C 

Polynomial 1 

0.1754 

0.2181 

0.04681 

Degree 2 

0.1370 

0.1343 

0.04677 

3 

0.1358 

0.1327 

0.03832 

4 

0.1323 

0.1207 

0.03668 

5 

0,1316 

0.1097 

0.035% 

Convective 

0.1371 

0.1530 

0.04683 

Shell 

0.1516 

0.1344 

0.04538 

Cloud Layer 

0.1444 

0.1347 

0,04556 

Intermediate 



0.05466 


Table III. Sunirary of Fitted Results 


For large N-M, the x 2 statistic, using (35), defined by (2x 2 )^ - 

(2N - - 1)\ is approximately Gaussian distributed with vanishing mean 

2 2 

and unit variance. Since we do not know o , we cannot enploy the x test 

2 

directly. However, the ai nature of the x distribution assures 

us that relative differences between model residual variances of only a 
few percent can be significant. 

In Case A, the hot areas of the North and South Eouatorial belts, 
the quadratic power series and the convective equilibrium model provide 
almost equally reliable fits, while all other models are much less probable. 
The convective equilibrium fit (19) gave the parameters b and c values of 
2.01 and 2.04 respectively. The value of b could be varied over a wide 
range (while that of c was adjusted in order to minimi re U for a given b) . 
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However, the value of c did not change significantly. Since c =2, the 
correspondence between the goodness of the quadratic and the convective 
equilibriun model fits is not unexpected. We discuss the physical dupli- 
cations of this result in the next section. 

The South Tenperate Belt is sure problematic . The Investigation 
of power series expansions for Case B reveals a preference for a fit of 
high degree, indicating significant structure. The convective equilib- 
rium model is clearly rejected. However, the other three-parameter models 
(the quadratic polynomial, the thin shell model and the cloud layer models) 
are equally likely, statistically! We can only conclude that the physical 
mechanism responsible for the behavior of the South Tenperate Belt is an 
amalgam of several physical processes or a region of transition between 
two physical processes . 

Finally, the intermediate flux region of the North Equatorial Belt 
allows for several models as possible mechanisms. The preferred degree of 
a polynomial fit is unity, as we have discussed earlier, while the addition 
of a quadratic term does not significantly change the results. The con- 
vective equilibriun model is viable and the associated parameter b can vary 
from 1.4 to 1.6 (while c =0.9) without significantly affecting the residual 
variance. Note that the Eddington approximation corresponds to b - 1.5 and 
c - 1.0. This is highly suggestive of flux conservation and radiative 
tlibriun. We also observe that the intermediate zone model of Terri le 
and Westphal is clearly rejected, while the thin shell and cloud layer 
models are approximately equally probable. The latter suggests that there 
is more thermal structure present than we normally associate with a state 
of radiative equilibriun. 
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Finally, let us consider how Lewis' (1969a , 1969b) model is 
consistent with these results. The hot spots have a brightness 
temperature of 250-255°K. This corresponds directly to the transi- 
tion \‘cg:Lon from aqueous acrncnia to ice, suggesting that a phase 
change coupled through the high opacity of aqueous armrnia (and the 
dynamic mixing that might arise fran this low-lying cloud layer) to 
corrective equilibrium may be present. The intermediate flux zone of 
the North Equatorial Belt is cooler and at a higher altitude. Since 
no phase transition is predicted and higher-level clouds are less 
likely to be dynamically coupled to what lies below, radiative 
equilibrium might be a reasonable approximation to the mechanism present. 
Finally, the South Tenperate Belt remains an enigma. Since it appears 
to be fairly hot, it could be intermediate in structure between the other 
regions. 


IV. CHEMISTRY CF PHASE TRANSITIONS 
In Section III, we argued for a phase transition in the hot spots 
of Jupiter’s equatorial belt. From the observation that the value of the 
convective equilibrium parsmeter c - 2.04 and the identification, from 
equations (18) and (19), that 

c - n (>~l)/u> , (47) 

(cf. Pollack a c* Sagan, 1965}, we can now estimate n, the term defined by 
(17) that shows the departure from a dry adiabat. We adopt n ’ 11.6. 

We now assume that the atmosphere of Juniter is 88 . 6 % hydrogen and 
11.2% helium, from Vfeidenschil ling and Lewis (1973). This is consistent 
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with solar abundances and with the 6 Sco occultation data (Elliot 
et al . , 1974) , Pioneer 10 ultraviolet photometer data (Carlson and 
Judge, 1974), Pioneer 10 infrared radiometer data (Orton, 1975) and 
the infrared spectnm determination by Houck et al. (1975) . A 
variation of 57, in our assuned He abundance will affect our results, 
through y, only about 17„. A simple calculation then shows that y - 1.42. 
Confining these results, we find that u - 1.68. The latent heats of 
water and anmonia are, respectively, 333.6 and 351 Joules per gram 
( International Critical Tables , 1928) , Now, neither water nor annonia 
will freeze at 250°K and the corresponding (Sagan and Salpeter, 1976) 

P - 2 bars (water freezes at a higher temperature and amaonia at a 
lower one). If they could, equation (17) would provide for mixing ratios 
of 0.0452 and 0.0494. These values are 5.35 and 42.8 times the estimates 
given by Weidenschilling and Lewis for solar abundance. Since, for 
clouds, we expect to find an excess of the spectroscopically active ma- 
terials, these excess values are not excluded. 

A mixture of water and aniiDnia, however, nroviJes for a large 
range of freezing points. (See Zemansky, 1968 and 0. stellan, 1971 for 
a discussion of eutectic curves and freezing mixtures.) The eutectic 
properties of aqueous ammonia solutions were investigated over half a 
century ago by Potsma (1920) and Elliott (1924) . At low temperatures , 
water does not readily dissociate in the presence of anmonia to form 
amroniira hydroxide. Anmonia can, however, form two hydrates, 
and NH^‘2H20 by hydrogen bonding. The two hydrates also exisL on the 
freezing point diagram. Therefore, depending on the strength of the 
initial aqueous anmonia solution, the sequence in which freezing takes 
place can be very conn lex. Moreover, unlike the laboratory situation, 
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the behavior in the Jovian atmosphere is considerably complicated by 
precipitation. Water is dense: than ice which is denser than frozen 
amim ia vMch is denser than liquid ammonia (Kuiper , 1952) . So , once 
the tenperature drops to that on the eutectic curve, one of the four 
active constituents 0 ^ 0 , NH-j^O, NH 3 ) will begin to freeze 

and either rise or sink faster than the mixture. The depletion of 
this constituent from the solution changes its concentration and lowers 
its associated freezing point. As the solution is buoyed higher by 
convection, it further cools mid loses more of one of its constituents. 

As a result, the freezing point can be smeared out over as nuch as 100°K 
(for the NH^ — system). Since the latent heats of fusion for 
water and ammonia are quite similar, and the bonding associated with the 
hydrates of anmonia are quite weak, the latent heats of the two hydrates 
should not be significantly changed and this picture remains unaltered. 

Water and anmonia vapors are relatively poor absorbers from 4.6 to 
5.1 vsn. Moreover, as we expect both anmonia and frozen anmonia to form 
above the water or ice clouds (from the above buoyancy arguments and 
Lewis, 1969b) , we expect that it would be very difficult to see down to 
the water or ice clouds at this wavelength. I Perhaps most water vapor 
that exists above the water clouds readily dissolves in the anmonia clouds 
and immediately freezes out. Since we expect large-scale moist convection 
to occur below this level (Gierasch, 1976) , the amount of water vapor 
present at lower levels will depend on Aether we are seeing a convective 
updraft or downdraft. In the case of a downdraft, the large-scale con- 
vective model predicts the presence of very little water vapor. The 
reaction in the expected absorption from water vapor in the downdraft 



would result in observations of much deeper and hotter levels In the 
Jovian atmosphere, The ecwbinatien, then, of a freeze-out mechanism 
at higher levels and convective downdrafts below could explain Larsen 
•t al. , 'a (1975) unexpectedly low water vapor abundance and high observed 
brightness temperature,} 

The possibility of hydrogen sulfide playir^ a major role in con- 
vective equilibrium models eamot definitely be excluded. Unlike water, 
R>S readily dissociates in ammonia. Moreover, all attempts to freeze 
such mixtures in laboratories have produced many expounds of snmonia 
and hydrogen sulfide. The eutectic curve for the NHj— H ? S system is in- 
completely la>ewn and appears to have at least one ireongruent melting 
point (Scheflan and McGrosky, 1932) ; in addition the relevant latent 
heats f«ve not been tabulated. For a reasonable estimate of the latent 
heat, the required mixing ratios are far in excess of that predicted from 
solar abundances. If hydrogen sulfide, however, is not well -mixed in the 
atmosphere, its role in the chemistry and coloration of the hot snots 
cannot be disccnnted (cf. Khare and Sagan, 1975). 

Although carbon monoxide has recently been detected in the 5 \m 
band in the Jovian atmosphere (Beer, 1975), it is tnlikely that 00 could 
be responsible for Terrile and Vtestphal ’ s observations. It is believes.! 
to be formed deep in the atmosphere (Larson et al , , 1978) . In the 
temperature range of interest, it can be expected to react with molecular 
hvvtrogen and to form methane (which is transplant at 5 un) and water . 
bbrtwvc , to fit the convective equilibrium wef-auiahat model , carbon 
monox tde would have to undergo a phase change near 250 C> K and would 
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necessarily be many tens or hundreds of times more abundant than solar 
values would surest. 

The analysis techniques of the present paper should be applicable 
to high-spatial resolution multi-frequency limb-darkening scans of 
Jupiter with the Voyager spacecraft (Hanel , et al. , 1978) — which can 
potentially clarify much about the lateral and vertical structure, 
chemistry and cloud constituents of the Jovian atmosphere. 



Vte are grateful for stimulating discussions to Peter Gierasch, 
E. E. Salpeter, William Forrest , Adam Devir, B. N. Khare, R. Treffers 
and Freeman Dyson. 

This research was supported by NASA Grant IJTR 33-010-082. 



29 


REFERENDA 

Abramowitz, M, and Stegun, I, A, (1965). Handbook of M athema tical 
Functions , Super intoickait of Dociments, Vfeshington. 

Aka Ike, H. (1969). Fitting autoregressive models for prediction. 

Am . hist . Statist . Math . 21, 243-247, 

Beer* R, (1975). Detection of Carbon Monoxide in Jupiter, Ajs, J. 

Lett . 200 , U67-L169. 

Carlson, R. W, mid Judge, D, L. (1974). Pioneer 10 ultraviolet photo- 
meter observations at Jupiter encounter, J. Ceophys . Res. 79 , 
3623-3633. 

Castellan, G. W. (1971). Physical Chemistry , Addison-Mesley, Reading, 
Mass. 

Chandrasekhar, S. (1960), Radiative Transfer , Dover, New York. 

Conrath, B. J. (1972). Vertical Resolution of Temperature Profiles 

Obtained from Remote Radiation Measurements , J . Atmos . Sci , , 29, 
1262-1271. 

Elliot, J., Vlassemrm, L. H. , Veverka, J., Sagan, C. , and Liller, W. 
(1974). The occultation of Beta Scorpii by Jupiter. 11. 'Hie 
hydrogen -hoi hm abundance in the Jovian atmosphere, Ast rop hysics 
J. 190, 719-729. 

Elliott, L. D. (1924). The Freezing Point Curve of the System Via ter - 
Airmen! a , J. Phys . Chem . , 28, 887-888, 

Forsythe, G. E. , Malcolm, M. E. and Holer, C. B. (1977). Computer 

Methods for Ma thema t ical Com p utation , Prent ice-Ha 1 1 , Engl evood- 
Cliffs, pp. 192-239. 

Foster, M. (1961). .An Application of the Wiener -Kolmogorov Smoothing 

Theory to Mat rix Inversion, J. Soc. Tndust. Appl . Math. 9, 387- 392. 



Gierasch, P. J, (1976). Jovian Meteorology; large -Scale Moist Compel ion , 
Icarus 29 , 445-454 . 

Gillett, F, C,, Low, F. J, and Stein* W. A. (1969). The 2.8—14 Micron 
Spectrum of Jupiter, Astrophy s . J. 152* 925-934. 

Goody, R. M, (1965). The structure of the Venus cloud veil, J. Geophys . 
Res. 70, 5471-5481. 

Graybill, F. A. (1976). Theory and Application of th e L inea r M odel , 
Duscbury Press, North Scituate, Mass. 

Hanel, R. , Gonrath, B,, Gautier, D. , Gierasch, P. , Kumar, S,, Kunde, V., 
lounvm , P. , Maguire, W. , Pearl, J., Pirraglia, J,, Pcrinamperum , C. , 
and Samuelson, R. (1978). The Voyager Infrared Spectroscopy and 
Radicrnetry Investigation, Space Sci ence Rev s . , in press. 

Hess, S. L. (1959). Introdu ction to The oreti cal Meteorol ogy , Holt, 

Rinehart and Winston, New York, 

Houck, J, R. , Sehaaek, D., Reed, R. A., Pollack, J. and Sunr.iers, A. 

(1975). Ji^iter: its infrared spectrum from 16 to 40 microns, 

Science 189, 720-722. 

T nt emaM onal Cri 1 1 cal Tabl es (1928) , MeGniw-Hi 1 1 , New York , 

Irvine, W. M. and Pol lack . J. B, (1968). Infrared Optical Properties of 
Water and Ice Spheres , Icarus 8, 324-360. 

Jenkins, G, M. and Watts, D. G. (1968) . Spectral Analysis and Its App lica- 
tions, Holden-Day, Son Francisco. 

Keay, C. S. L. , Low, F. J., Rieke, G. H. and Minton, R. B. (1973). High- 
Resolution Mips of Jupiter at Five Microns, Astrophys. J. 183, 
1063-1073. 



31 


Khare , B, N. and Sagan, C, (1975), Cyclic Qctatomie Sulfur; A Possible 
Infrared and Visible Gtromophore in the Clouds of Jupiter, Science 
189, 722-723, 

Kuiper , G, P, (1952), Planetary Atmospheres and Their Origin, in The 
Atrospheres of the Earth and Planets (Revised Edition), G. P, 

Ruiper, editor, U. of Chicago Press, Chicago, 

Larson , H, P,, Fink, U., Treffers, R, and Gautier, T, N, III (1975). 

Detection of Water Vapor on Jupiter, Ast ronhy s , J. 197, U 37 -LIAO. 

Larson, H, P,, Fink, U. and Treffers, R, R, (1978). Evidence for CO 
in Jupiter’s Atmosphere fron Airborne Spectroscopic Observations 
at 5 Microns, J. 219, 1084-1092. 

Lawson, C. W, and Hansen, R. J, (1974). Solving Least Squ a res Pr o blems , 
Pm\tice-Hall, Ehglewood-Gliffs, 

Lewis, J. S. (1969a), Observability of Spectroscopically Active Ccnpounds 
in the Atmosphere of Jupiter, I caru s 10, 393-409. 

Lewis, J, S, (1969b), Tlx? Clouds of Jupiter and the and NH^-l^S 

Systems, Tearus 10, 365-378. 

Luenberger, D. G. (1973), Introduction to Linear and Non-Linear Pro- 
graming, Addison-Woslev, Reading, Mass, 

Ntwivrn, W. I. (1975). Infrared Lircb Darkening', of the Venus Atmosphere, 
Ica rus 26, 451-456, 

Orton, G, S, (1975). The thermal structure of Juoiter. I., Icarus 26, 
125-141. 

Oi'UTn , G. S. (1977). Recovery of the Moan Jovian Temperature Structure 
from Inversion of Spectrally Resolved Thermal Radiance Data, 

Icarus 32, 41-57. 



32 


Parker, R. L, (1977). Understanding Inverse Theory, Am. Rev. Earth 
P lane t, Set . , 5 , 35-64 . 

Pollack, J. B, and Sagan, C. (1965). The Infrared hirb -Harkening of 
Venus, J. Geophvs. Res. 70, 4403-4426. 

Potsma, S. (1920). Le Systems Anvnoniaque -Eau, Ree. Tray. Claim. , 39 , 
515-536. 

Ralston, A. (1965). A First Coiyse in llumeriral .Analysis, McGraw-Hill, 
New York. 

Robertson, C. W. , Downing, H. D. , Camutte, B. and Williams, D. 1,1975). 
Optical constants of solid antnania in the infrared, J. Opt. Soc, 
Am. 65, 432-7*35. 

Robertson, C. W. and Williams, D. (1971). Inmbert Absorption Coeffi- 
cient of Water in the Tnfrared, J. Opt, Soc. Am. 61, 1316-1320. 

Robertson, C. W, and Williams, D. (1973). Optical constants of liquid 
amwnia in the infrared, J. Opt. Soc. Am. 63, 188-193. 

Sagan, C. (1971). 'Dae Solar System Beyond Mars, Space Science Revs. 11, 
73-112. 

Sagan, C. and Sal peter , E. E. (19/6). Particles, Enviroi intent s and llv- 
not helical Ecologies in the Jovian Atmosphere, Astrophvs. J. 

Suppl . 32, 737-/53. 

Sagan, C. , Veverka, J,, Washerman, i. , Elliot, J, , and biller . W. (197*4). 
Jovian Atmosphere: Structure .and Corposition lx j twen the Turbo- 

pause and the Mesopause, Science 187*. 901-903. 

Seheflan, b. and MeCrosky, C. R. (1^33), A Phase Snub of the System 
Anisonia -Hydrogen Sulfide, J. .Auer . (hem. Soc. 54, 193-202. 

Terri lc, R. J. ami Weslphal , J. A. (l l V7). The Vertical Cloud Structure 
of Jupiter from 3 tmi Measurement s , Te'arus 30, 2/4-281. 


Weiden schilling, S. J. and Lewis, J. S. (1973) Atmospheric and Cloud 
Structures of the Jovian Planets, Icarus 20, 465-476. 

VJestphal, J. A. (1969). Observations of localized 5-Micron Radiation 
from Jupiter, Astrophys . J. 157, L63-L64. 

Westphal, J. A, , Matthews, D. and Terrile, R. J. (1974). Five-Micron 
Pictures of Jupiter, Astrophys . J. 188 , L111-L112. 

Zemansky, M. W. (1968). Hea t and Thermodynamics , McGraw Hill, New York. 



33 


APPENDIX I 


Orton (1977) enployed Conrath's (1972) formulation of the Backus- 
Gilbert generalized inverse theory. In the discrete case in the absence 
of noise, the method can be stated quite succinctly. Consider the 
problem of best approximating some function AT(x) given m observations 
AIj defined by 

*t 

At, - / k. (x) AT(x)dx, I - 1, ... ,m . (1-1) 

0 

(We adhere strictly to Conrath's notation. The discussion that follows, 
however, is independent of the choice of the limits of integration, pro- 
vided they are finite.) We then wish to construct an aoproximate inverse 
ATO^ from a linear ccnfcination of the observed data, namely 

AT(xJ - ? a. (x) AI. (1-2) 

i-1 1 1 


where we have yet to specify a rule for selecting a^(x). If we define a 
function A (x,x') by 


A(x,x') * I a. (x) k, (x') , (1-3) 

i-1 1 1 

we see that 

AT(x5 - rACx.x'jATCx’)*' . (1-4) 

0 


The function A(x,x') , ideally, should tend to a Dirac 6-function. In 
practice, it will have a finite width or spread and tends to smooth AT(x). 
For this reason, it is called an "averaging kernel." In order to estimate 
the width of A(x,x'), we define a "spread finction" s(x) by 
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•00 - 12 r Ot-x') 2 A 2 (x,x'> dx* . (1*5) 

0 

We then perform a variation on s(x) in order to minimize the ’’spread,” 
subject to the normalization constraint 

*t 

l-rA(x,x')dx’ . (1-6) 

0 

The factor 12 in equation (1-5) is Introduced so that, if A(x,x') is a 
rectangle of unit area and width w, s (x) ■ w. This variation can most 
sinply be achieved using Lagrange nultipliers (Conrath , 1972) . 

The method is conceptually attractive since it will provide the 
estimate of AT(x) with vrtiat seems to be the best possible resolution. 

Although the method can be of significant value in certain applications, 
it suffers from important raathanatical shortcomings which can seriously 
affect its performance . In the absence of any further information about 
the physical processes involved, two mathematical principles must be em- 
ployed when devising an ad hoc inversion schane. First, all available 
information rust be incorporated into the method so that the solution ob- 
tained reproduces the available data. Second, the method must yield in- 
creased resolution over other approaches. 

The variational procedure employed in estimating a^(x) never uses 
the information obtained by the observations (1-1) as equations of constraint. 
As a result, the approximate inverse AT(x5 will not, in general, reproduce 
the observed data AI^. By not fully introducing the observational informa- 
tion available into the method, we compound our ignorance of the solution. 
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The second problem arises In the determination of the resolution 
or, alternatively, the spread, Although, s(x) reproduces reasonably 
well the width of a ranker of functions , it can give spurious results 
when applied to some degenerate kernels. Consider, for exanple, a 
hypothetical experiment there we measure AI^ define i by 


1 

A If - / P^CxMTCxJdx, i - 1 , . . . ,m (1-7) 

there P^(x) is the i^ Legendre polynomial. Let us select a^(x) and, 
therefore, “ATOcJ’ to be given by 

a i^ x > “ ^ P i-l (x) 


( 1 - 8 ) 


at(xT 


m 21-1 

E ■ • ’y — P^( x)AI^ 


(In many circumstances , we customarily make this choice of expansion 
since it provides the best approximation, in an integrated least-squares 
sense, to a given function.) Although Conrath's variational procedure 
would not make this identification for a^(x), it is instructive to consider 
the spread function that results from this choice. 

Now, the averaging kernel becomes 


A(x.x’) 


E „(x) P.(x') 

1*0 L 1 * 


(1-9) 


(We know, incidentally, from the completeness relation for Legendre 
polynomials that this kernel "tends" to a Dirac ^ -fine t ion. ) Using the 
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recurrence relation 


<*+l)P l+1 (x) - (21+DxP^Cx) + iP^OO - 0 


<I-10) 


we obtain the Christoffel-Darboux identity 

Then, using (1-5) with integration limits *1 to 1, we observe that 

r- n 


(i-ll) 


S(x) - 6m 2 


D *■ 

. m-l(x) . _ 

L“£*r + “5T"J 


p m 2 < x ) 


(1-12) 


and the average spread <s(x)> is 


<s(x)> * hf s(x) dx 
-1 

12m? 

- — 

4mM 


(1-13) 


This result shews that the average spread increases as we add more terms 
(and corresponding data points) and, as m approaches infinity, reaches a 
limiting value of 3 (vMch is larger than the region over which we are 
calculating the spread) . The measure of spread that we use nust shew in- 
creased resolution as we increase the number of data points and, in the 
limit of an infinite amount of available information or data, nust tend to 
zero. The spread function of (1-5) is inconpatible with these conceptual 
requirements . 
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Although we have shown that the spread f ^mction (1-5) can be a 
misleading Indicator of resolution, it is irrortant to understand in 
practical terms why this is so. The morphology of a typical averaging 
kernel is characterized by a central peak, for x near x’ , and sane kind 
of “sicmobe” structure. Theoretical kernels (e.g. , rectangles, Gaussians, 
etc.) for which (1-5) is a reasonable measure o'* spread have no sidelobes . 
The kind of kernel more likely to be encountered in practice will have a 
complex sidelobe structure. If these sidelobes do not decay nuch faster 
than (x-xT 1 , they will provide a significant if not dominant contribu- 
tion to the spread function (1-5). In practice, the highly oscillatory 
character of the sidelobes tends to cause can: dilations aid diminish any 
sizeable contribution. However , since (1-5) ccnta^is the square of t> 
kernel, this cancelling feature of the sidelobes is lost and their effects 
are grossly exaggerated. 

This form of generalized inverse theory, then, has two serious draw- 
backs. It will nor. reproduce the given data and can provide a very spurious 
estimate of the resolution of the result. For completeness , we cite an ap- 
proximate inversion fornila (see, for example, Foster, 1961) that satisfies 



that is, ary error in our approximate inverse cannot be seen from available 
observational data. We define a matrix C by its i, j comnonents 


C. . - / t k. (x) k. (x) dx i , j ■ 1, . . . ,m 

i.J 0 J 


(1-15) 
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Then, 


a* 00 - E 
1 j 


i,J k j°° 


A(x,x’) - ^ kjCx) k. <x') 


(1-16) 


m 

T(x) « l 

t.J-1 


,»1 


l,j X i k j (X) 


wtere |C _1 1 . . denotes the i, j tk conponent of the Inverse matrix to C. 
J 



